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ABSTRACT 

We discuss the extent to which photometric measurements alone can be used to identify Type la 
supernovae (SNIa) and to determine the redshift and other parameters of interest for cosmological 
studies. We fit the hght curve data of the type expected from a survey such as the one planned 
with the Large Synoptic Survey Telescope (LSST) and also to remove the contamination from the 
core-collapse supernovae to SNIa samples. We generate 1000 SNIa mock flux data for each of the 
LSST filters based on existing design parameters, then use a Markov Chain Monte-Carlo (MCMC) 
analysis to fit for the redshift, apparent magnitude, stretch factor and the phase of the SNIa. We find 
that the model fitting works adequately well when the true SNe redshift is below 0.5, while at z < 0.2 
the accuracy of the photometric data is almost comparable with spectroscopic measurements of the 
same sample. We discuss the contamination of Type Ib/c (SNIb/c) and Type II supernova (SNII) on 
the SNIa data set. We find it is easy to distinguish the SNII through the large mismatch when 
fitting to photometric data with la light curves. This is not the case for SNIb/c. We implement a 
statistical method based on the Bayesian estimation in order to statistically reduce the contamination 
from SNIb/c for cosmological parameter measurements from the whole SNe sample. The proposed 
statistical method also evaluate the fraction of the SNIa in the total SNe data set, which provides a 
valuable guide to establish the degree of contamination. 

Subject headings: cosmology: theory — distance scale — large-scale structure — supernovae: general 



1. INTRODUCTION 

The cosmological applications of luminosity-distance 
measurements to Type la supern o vae (SNela) are 
now well known (iRiess et al.l 119981 : iPerlmutter et al.l 
119991 : iLeibundgutl 1200 1[ ). While the current sample 
of SNela-b ased distances ar e limited to a few hun- 
dred SNe fAstier e t al.' '20 0I IWood-Vasey et all [20071: 
Kowalski et all |2008: Hick en et al.l 120091: iKessler et al.l 
2009^, future surveys are now planned to increase the 
sample size to a few thousand or more that could poten- 
tially allow a few percent accurate dark energy equation 
of state measurements in several redshift bins between 
< z < 1 (see e.g. iHowell et all 120091: iSarkar et all 
[2008h . The main challenge for constructing large samples 
are likely to be spectroscopic follow-up measurements to 
identify if each supernova detected in a photometric mon- 
itoring campaign is Type la and to establish the redshift 
of that supernova. 

In addition to the planned space-based programs such 
as the Joint Dark Energy Mission (JDEM) ^, in the near 
future, there will also be several ground based photo- 
metric surveys for cosmological measurements and other 
astronomical studies. These include the Dark Energy 
Survey (DES)^, the Pan-Starrs survey ^ , and ultimately 
the Large Synoptic Survey Telescope (LSST) which 
plans to monitor a large area of the sky every few days 
leading to a large sample of transient sources includ- 

^ http://jdem.gsfc.nasa.gov/ 
^ http://wwvir.darkenergysurvey.org/ 
^ http:/ /pan-Starrs. ifa.hawaii.edu 
* http://www.lsst.org/ 



ing supernovae. Given the large size of the samples of 
SNe expected, it is highly unlikely to have spectroscopic 
follow-ups for all or even a large fraction of them. Due 
to this limitation it appears challenging to obtain cos- 
mological measurements with the SNIa seen by LSST. 
Since LSST is likely to detect a few hundred thousand 
or more SNe per year, it would be highly desirable to 
identify whether a given SN as Type la or not, and to 
extract useful parameters such as redshift and luminos- 
ity with photometric data alone. If reliable techniques 
could be established, then even with a large degradation 
in accuracy for individual data compared with the case 
where spectroscopic data are also available, given the 
large number statistics expected, one could still aim to 
achieve a good measurement of cosmological parameters. 

In this spirit we pursue a study to establish the extent 
to which photometric data from a survey like LSST can 
be used to identify SNela and to measure the cosmolog- 
ical parameters. We do this by fitting the photometric 
light curve data with sampling and errors consistent with 
LSST. Our mock SNe samples also include core-collapse 
supernovae and we vary the fractions expected based on 
the current rate estimates of various types of SNe. Our 
MCMC analysis are focused on a joint parameter estima- 
tion including the redshift, apparent magnitude, stretch 
factor and the phase of the SNIa. We find that the model 
fitting works adequately well when the true SNe redshift 
is below 0.5. At 2; < 0.2, photometric data of the type 
expected with LSST provide an accuracy that is close to 
the case when spectroscopic measurements are also avail- 
able, with the redshift determined separately from spec- 
troscopic data leading to one less parameter in MCMC 
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fits than photometric light curves. 

We also focus on the contamination of Type Ib/c 
(SNIb/c) and Type II supernova (SNII) on the SNIa data 
set. We find it is easy to distinguish the SNII from SNIa's 
through the large mismatch in the fitting to la light 
curves. This is not the case for SNIb/c and they pro- 
vide the main contamination to la measurements. In 
addition to a cut in values, we implement a statistical 
Bayesian estimation method to reduce the contamination 
from SNIb/c in the subset of SNe sample selected for cos- 
mological measurements. This technique also establishes 
statistically the fraction of the SNIa in the total SNe 
data set, which provides a valuable guide to the degree 
of contamination from Ib/c's. 

We employ the filter functions as currently publicized 
by the LSST team in addition to survey parameters 
outlined in Ivezic et al. (2008). We note that while 
our work is focused towards a survey like LSST, oth- 
ers have also consider the use o f photometric data alone 
for SNe distance me asurements f Johnson & Crotts'l2006|; 
Sullivan et aLll200a Kuznetsova & ConnoUv 2007; Wang 
2007t IKim fc Miquell i2007l: iZentner fc Bhattachar^ 
20091 ^. 

The discussion is organized as follows. In the next 
Section we describe our procedure to simulate SNela data 
in a survey like LSST and move on to discuss our six 
parameter model fits to the multi- wavelength light curves 
from a large mock sample using a MCMC analysis. In 
Section 4, we discuss the contamination from Type II 
and Ib/c SNe to la photometric samples and a way to 
statistically reduce the contamination from Ib/c's using 
a technique that implements the Bayes theorem. 

2. SIMULATING SNIA OBSERVED FLUX DATA 

In this Section, we describe the process to generate the 
various SNe data. We first discuss the observed SNIa 
mock fiux data. 

2.1. The Mock Light Curve 

The apparent observed fiux from a supernova at z can 
be written as the convolution of the spectral energy dis- 
tribution (SED) and the transmission function of the 
telescope, 

fobs = J Tx{Xobs)SE^{^obs,tobs, S, Eg°^Y, Eg^)dXobs, 

where Tx{Xobs) is the filter response for band X, Xobs is 
the observed wavelength, tobs is the observation date, s 
is the stretch factor, and Eg°fy and E^]^ are the color 
excess for the host galaxy and the Milky Way respec- 
tively. 

The transmission functions used in our analysis for 
S-ban ds of LSST are shown in Fig[T] (|LSST filter draftl 
[2005h .^ We also plot rest-frame SEDs for Type la SNe 
at z = 0.2, 0.5, and 0.8. These SED templates are from 
Nugent et al. (2002) and they cover the spectral wave- 
length from 1000 to 25000 A in rest-frame days from —20 
to 70 with respect to the B-band maximum light day. 
The SNe fiux for the epochs before -20 are set to be zero. 

^ We note that there are several filter designs for LSST including 
a scenario involving 6 filters. Here, we focus on the 5-band case 
with simple filters. 




Fig. 1. — The LSST filters used in our analysis. From left to 
right are g, r, i, z and y-band functions. The Type la SEDs of SNe 
for 2 = 0.2, 2 = 0.5 and z = 0.8 at the day of rest-frame B band 
peak magnitude are also shown. The flux is in an arbitrary unit. 

There are now several techniques to parameterize the 
SNIa light curves, such as the 1 5-da y decline after the 
B-band maximum light Atoi5 (jPhillip^ |1993) and the 
multicolor light curve shap e method (MLCS and the 
upda te version MLCS2k2) (IRiess et all 11991 IJha et all 
[200l . In this paper, we calibrate the SNIa light curve 
with the time-scale an d stretc h factor relation following 
the works of Perlmutt er et all ^997, 1999). By stretch- 
ing and compressing the time axis around the rest frame 
B-band maximum light day, this method can fit the ob- 
served light curve very w ell using the light curve template 
(|Goldhaber et al.ll200lD . Then this SED can be re-scaled 
by the apparent, unextincted B-band peak magnitude 

ms = Mb + 51ogiorfL(z, 0) + 25 ~ a{s ~ 1) + Am, (2) 

where Mb is the B-band absolute peak magnitude, cIl 
is the luminosity distance which is a function of the red- 
shift z and a broad set of cosmological parameters de- 
noted by 9 and a is the coefficient of the relation be- 
tween s and Tor. Here we ta ke Mp — —19.3, a = 1.5 
(|Knop et al.ll2003l : IAstier et al.ll2006h . and the set of cos- 
mological param eters with flrnp = 0-27, ^Iao = 0.73 
and Hq — 0.71 (jKomatsu et al.l [2009f) where f)mO and 
r^Ao have the usual meaning with the present-day matter 
and dark energy density parameters and Hq is the dimen- 
sionless Hubble constant. Besides, we also consider the 
dispersion Am of the rest-frame B-band peak magnitude 
after the calibration of the stretch factor. The B-band 
filter wc use is from the Johnson- Morgan system (Bessell 
[1990, 2005). Also, the time scale of the SED is calibrated 
by the stretch factor s, which is assumed to be available 
from —15 to 35 around the B-band maximum luminosity 
day (jAstier et al.l 120061 ). We note that there is also an 
intrinsic color scatter, crf^^ (standard deviation), that 
should be taken into account when producing mock light 
curves. Based on prior work, we find this uncertainty to 
be small with a value of ~ 0.05 mag (Ph ilhps et al.l[T99l 
iJha et aLll2007f) . Hence it would not affect our results 
much, so that for simplicity we don't consider it here. 
During the transit, the supernovae light will be partly 
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absorbed by the dust of the host galaxy. We employ 
the reddening law of Cardelli et al. (1989) with Ry — 
3.1, from infrared to far-ultraviolet (0.3/im~^ < a; < 
10/im~^, where x = 1/A). For the optical to near ultra- 
violet wavelength range (l.l/Ltm~^ < x < 3.3^m~^), we 
use an updated version for extinction given by O'Donnell 
(1994). The latter uses the same analytical form for ex- 
tinction as Cardelli et al. (1989) but with values of the 
fitting parameters revised slightly from the previous ver- 
sion. The level of extinction we assume here is consistent 
with the one measured recently by Menard et al. (2009) 
corresponding to large angular scales based on galaxy- 
QSO cross-correlation in SDSS. 

Since the SNe are at a different redshift than the ob- 
server, the spectrum is redshifted for both wavelength 
and the phase, i.e. A' = A(l -I- z) and t' = t{l + z). We 
also a pply an extinction associa.ted with dust in the Milky 
Way (jBurstein fc HeiIei[T98l FSchleeel et al.lll998r) . and 
assume that we have a perfect measurement of Eg]^. 
This assumption has no effect on our final conclusions. 
While the extinction of the Milky Way have different val- 
ues for different sky regions, we do not have any infor- 
mation on the exact field selection of future SNe surveys 
from ground. Thus, we do not account for sky variation 
of extinction and simplify by just taking an average value 
with « 0.03 and Ry = 3.1 for the extinction law. 

Finally, the spectrum is integrated with the LSST filters 
to get the mock light curve sampling in each of the 5 
LSST filters. 

Since the mean redshift of the SNe detections with 
LSST main survey is expected to be about 0.5, and the 
deeper, but smaller, survey can potentially detect SNe 
out to ~ 1, we choose the redshift range from 0.01 to 1.1 
when making mock SNe samples. 

When creating large samples, we assume the flat 
ACDM model with Qmo = 0.27 and ho = 0.71. Then, 
the redshift z, stretch factor s, the extinction of the 
host galaxy Eg°fy and the magnitude dispersion Am 
are generated from the Gaussian distribution with trun- 
cated tails as follow: 0.01 < z < 1.1 with z = 0.5 and 
= 0.4, 0.6 < s < 1.4 with s = 1 and = 0.3, 
-0.1 < E'^°f\, < 0.3 with E'^% = 0.0 and as = 0.2 

and -0.3 < Am < 0.3 with = 0.0 and aAm = 0.17. 
This extra disp ersion acts as an extra source of noise in 
our mock data (Sullivan et al.' 2006: Har nuv et all 119951 
[1996; Phillips ct al. 1999; Guy ct al. 200^. 

2.2. The Photometric Error and The Cadence 

The photometric error we use for LSST comes from 
Ivezic et al. (2008) and takes the form of 

^phot ^ sys ^ zero ^rand^ } 

where Usys is the systematic photometric error which is 
designed to be very small (< 0.005mag). azero is the 
absolute photometric error that we set to be Ozero = 
0.02mag (|Astier et all l2006t iSullivan et all [20061 ). We 
note that, in practice, there is only one zero-point real- 
ization in any given experiment that is applied to all su- 
pernovae. This would result in a n on-diagonal covarianc e 
matrix for the distance modulus (|Kim fc Miauel|[2"006[ ). 
However, since the inclusion or non-inclusion of this co- 
variance does not change the principles of our method- 



ology, for simplicity, we ignore this correlation here. We 
do suggest that it must be considered in an analysis of 
real data. 

In equation (3), cr^(j„^ is the random photometric error 
for point sources given by 

CTl,^, = {Qm-^)x + ^x\ (4) 

Here 7 is a parameter related to the sky brightness and 
readout noise, among others, and x = ioO '*('"~"'s), 
where m is the magnitude and ms is ha depth for a de- 
tection of a point source in each of LSST bands. The ms 
is a function of the sky brightness, the seeing, the expo- 
sure time, atmospheric extinction, the airmass and the 
overall throughput of the instrument. All of the value of 
these parameter can be found in Table 2 of Ivezic et al. 
(2008). 

We randomly generate the first observational day from 
-20 to 35 rest-frame days to ensure that we always have 
enough data to establish the stretch factor. We next 
randomly select the data point to occ ur every 3 or 4 days 
based on the cadence of the LSST (jlvezic etall [20081. 
Finally, about 1000 mock SNIa flux data are generated 
for each of the five filters. 

In Fig(2J we show the examples of the mock light curves 
in g, r and i bands at different redshifts. The mock flux 
is created from the Gaussian distribution with the mean 
on the light curve. Here we set the first observe-day 
t?^'^ — —10 in the observer- frame. 
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Fig. 2. — The examples of LSST mock SNIa light curves and 
observational data. The solid blue, dashed cyan and dotted green 
lines are the g, r and i band light curves respectively, and the first 
observe-day is set at -10 day in the observer-frame. The flux is in 
an arbitrary unit. 



3. FITTING THE LIGHT CURVE 

There are six light curve parameters that we hope to 
extract from multi- wavelength light curve fitting. These 
parameters are the z, ms, s, E'^°_^y, Am and io*^"' (i.e. 
the rest- frame date for the first observe-day). The 
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statistical method is employed here with 



t bands 

EE 



£ohs 
Jij 



-obs 



where /°/^ ^ 



(5) 

fj', and (T°j^ = cjphot are the observed, theo- 
retical flux and observed error for the observe-day ti and 
band j, and Tj is the transmission of band j. The sum- 
mation goes through all bands and days with observed 
samplings of the light curves. 

3.1. The Markov Chain Monte Carlo Technique 

The best-fit value for each light curve parameter usu- 
ally can be found using the nonlinear least-squares fitting 
technique (e.g. sullivan et al. 2006). Here considering 
the number of the parameters, the efficiency and the ac- 
curacy, we would like to employ the MCMC technique 
to perform the fitting process. This method does not 
require to assume a Gaussian distribution for the like- 
lihood, and it is easy to perform the marginalization 
over other parameters when quoting error for one pa- 
rameter. Most importantly, it is very efficien t for the 
multi- pa r amete r fitting (Neil 199^ [Lewis fc Br idle 20021 

20041: IGong fc Chenl 



MacKavl [200l [P oran fc MucUct 
20n7tlTrotta 1120081 ). 



Our purpose is to estimate the posterior probability 
P{9\D) for the parameter set 9 given the observational 
data set D. Based on the Bayes theorem 



p{e\D) 



L{T>\e)P{e) 



(6) 



where C{T>\9) ^ e^^ is the likelihood which denotes 
the probability to get D given the parameters 9, P(9) is 
the prior probability for 9 and -P(D) is the normalization 
factor which would not affect our analysis here. 

The Metropolis-Hastings algorithm is applied in our 
MCMC technique to decide if a new point should be ac- 
cepted by an acceptance probability: 



.(9 Ifl Wn.inJ:^felM^!i^ 1 



£(D|g„+i)q(g„|g„+i) 
/:(D|0„)q(0„+i|0„) 



(7) 



(8) 



where q(0„+i|0„) is the proposal density to propose a 
new point 0„eu) given a current point 0„ in the chain. 
Here we assume uniform prior probabilities for the pa- 
rameters which is canceled in Eq.Q. If a = 1, the new 
point 9new is accepted; otherwise, the new point is ac- 
cepted with probability a. This process are repeated un- 
til a new point is accepted, and then we set 9n+i = dnew 
Also, we set a uniform Gaussian-distributed proposal 
density for every point, so that it is independent of the 
position on the chain, i.e. q{9n+i\0n) — ^i9n\9n+l), we 
then have 



(9) 



Since the proposal density determines the step size of the 
MCMC process, it is closely related to the convergence 



and mixing of the chain. Here w e adopt the adaptive ste p 
2 size Gaussian sampler given by (jPoran fc Muelleill2004D . 
The criterion of the con vergence we use was described in 
' iGelman fc RubinI (I1992D. and after convergence we freeze 
the step size (jDoran & Mueller 2004). 

The ranges of the parameters when we run the MCMC 
are set as follow: z e (0, 2), tob G (10, 30), s e (0.5, 1.5), 
Eb°-v G (-0.5,0.5), Am g (-0.5,0.5) and tl""^ € 
(-20,40). For each mock SNIa, we take about 10000 
chain points to illustrate the probability distribution of 
the parameters after the burn-in and thinning process. 

3.2. The Light Curve Fitting Results 

In Fig 131 we compare the input redshift of each of our 
1000 SNela in the simulation with the redshift obtained 
from MCMC fitting of SEDs to the multi-wavelength 
light curves. We find that when z < 0.2 the SNela light 
curves are adequately sampled with enough accuracy to 
allow good redshift estimates along with other parame- 
ters, with uncertainties as small as 0.001. Such an error 
is comparable with the spectroscopic measurements, and 
even if the spectroscopic measurements could provide 
a higher precision on the measure ment of the redshift, 
in an y case th e unknown bu lk flows (jCoorav fc Caldwell 
120061: [Zhang fc Chen |[200l) would produce an error on 
the redshift at this level. For the medium redshift 
0.2 < z < 0.5, the estimated redshift is still useful but 
the the uncertainty is about 0.1. For z > 0.5, the appar- 
ent magnitude becomes large, and since aphot ^ 10° *™, 
the redshift errors increase quickly with increasing red- 
shift and can reach ~ 1. The limitation at high redshift 
is also due to lack of near-IR photometric coverage and 
addition of IR bands beyond the z-band will improve 
photometric determinations when z > 0.5. 

The residuals and la errors for the total six fitting 
parameters in the MCMC analysis are shown in FigHl 
Similar to the redshift, the multi-wavelength light curve 
model fitting leads to parameter accuracies that are re- 
markably accurate when z < 0.2, except for Am as it 
acts as an extra source of noise independent of the red- 
shift. As shown in FiglH over the whole redshift range 
studied out to z of 1.1, the dispersion of the fitting z, 
niB, s, Eg°_fy, Am and t^'^'^^ are mainly less than ±0.3, 
±0.4, ±0.2, ±0.2, ±0.4 and ±4, respectively. 

3.3. The Constraints on Cosmology 

To establish the overall effect of the uncertainty from 
photometric redshift for cosmological studies, we also 
generate 1000 SNela with spectroscopic redshifts Zgpec 
with the LSST photometric error aphot, i-e. we just fix 
the redshift and only model fit the other five parameters. 
The Hubble diagram for the spectroscopic and photomet- 
ric cases are shown in FiglB] Only one sixth of the whole 
data are shown on each figure. 

We use the MCMC approach to fit the cosmological 
parameters from the two Hubble diagrams. Two cosmo- 
logical scenarios are considered, that first one is ACDM 
with non-fiat geometry and the second is wCDM with 
the time-evolved equation of state for the dark energy 
with w{z) = Wo + wiz/{\ + z). 

In Figl?! we show the contour maps of il„jo vs. JIao 
and wq vs. wi withlcr and 2a errors. As can be seen, 
the la contours using Zphot and its error in cosmological 
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Fig. 3. — The intrinsic redshift of each of the 1000 mock SNcIa 
hght curves compared to the photometric redshift estimated with 
multi-parameter MCMC fits to multi-wavelength light curves. The 
best fit value and la errors are shown. The redshift estimation 
is remarkably accurate at low redshifts when z < 0.2 with errors 
comparable to either the spectroscopic measurements of redshift or 
theoretical uncertainty in the redshift coming from peculiar veloc- 
ities and bulk flows, among others. At z > 0.5, the fitted redshift 
errors are significantly larger because of the large photometric er- 
rors CFphot- 
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Fig. 4. — The residuals for the six parameters in the MCMC 
analysis with Itr errors for the 1000 SNela in the mock sample. 
The results are pretty good for z < 0.2, except for Am since it can 
be seen as the noise and is independent on the redshift. 

parameter fits nearly overlap with the 2a contours of the 
case where redshift is known precisely usingZgpec- Also, 
we find little deviation for the directions of the main axis 
of the contours for the two cases. Thus, for a survey such 
as those planned for LSST, we effectively find a factor of 

2 degradation in parameter uncertainties when using 
the SNIa sample with only photometric redshifts com- 
pared with one with also spectroscopic redshifts. 

In terms of the dark energy figure of merit that involves 
the inverse area of the wq vs. wi ellipse, photometric SNe 
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Fig. 5. — The distribution of the fitting value minus the actual 
value for each light curve parameter. The number has been nor- 
malized. 




0.5 




'phot 



Fig. 6. — The Hubble diagram for the Zspec and Zphot simulations. 
The data points on each figure are just one sixth of the whole data 
sets. 

samples lead to a factor of 4 degradation compared to 
spectroscopic sample. This difference, however, is likely 
to be a minor issue: compared with the planned SNela 
surveys which will involve spectroscopic measurements of 
a few thousand SNe per year, photometric only surveys 
such as the one with LSST will produce a sample of a 
few hundred thousand SNe. Moreover, we have to note 
that for a real survey the sample is always magnitude- 
limited, so that some z > l.l objects could contaminate 
the 2: < 1.1 sample and lead to a bias in the fitting results. 
Given that we cannot quantify this bias fraction, we don't 
include this effect in our analysis. In an upcoming paper, 
we hope to implement a new technique to account for 
such biases in large SNe samples. 

4. REMOVING THE CONTAMINATION FROM SNIB/C 
AND SNII 
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Fig. 7. — The contour maps for QrnO vs. Qao (left) and uiq vs. 
uii (right). The la and 2cr errors are shown. The red solid and 
blue dotted contours are for photometric and spectroscopic rcdshift 
simulations respectively. 

In a pure photometric survey such as the one with 
LSST without spectroscopic measurements to identify if 
each of the SN is Type la or not, in addition to the 
error in the measurement of redshift, the photometric 
SNe samples would also be contaminated by core-collapse 
supernovae. Here we consider the contamination from 
Type Ib/c (SNIb/c) and Type II supernova (SNII). 

4.1. Estimating The Contamination 

To estimate the level of contamination, we create mock 
samples of light curve data for SN Ib/c and SN I I. Th e 
spectral templates for SNIb/c is from lLevan et al.l()2005D . 
and for SNII we use the temp lates of SNIIP an d SNIIL 
given bv lGilliland eFaLl ([19991) and lBaron efall ([2004). ^ 
Here, for simplicity, we consider primarily the SNIIP 
and SNIIL. The S NIIn which have "unusual" progenitors 
(|Mobberlevl [20071) may be an import ant contamination 
for the SNIa (jPoznanski et al.ll2007bl ). We may discuss 
these objects in future work. We set the percentage of 
SNIIP and SNIIL are 50% and 50% respectively for the 
SNII sample. The mock flux data of the SNIb/c and 
SNII samples are generated with the same procedure as 
those used to generate the SNIa mock data in § 2. 

Since the core-collapse supernovae are intrinsically 
fainter than the SNela and have no magnitude-phase 
relation, we take the absolute peak magnitude from 
[Richardson et all ()2OO20 and the Gaussian distribution 
with z = 0.4 for SNIb/c and SNII, and then set s = 1 
when mimicking their B-band peak magnitude. Also, 
given that the core-collapse supernovae are usually found 
in star forming regions, they are expected to suffer more 
extinction from the host galaxy. We set —0.2 < E'jffy < 

0.4 with E'^°_fl, = 0.0 and as = 0.3. Once simulated, we 
continue to use the SNIa SED light curves to fit them 
just as we did in § 3. 

The distribution of the difference of Xmm (i-^- relative 
X^) for the SNIa, SNIb/c and SNII are shown in FigB 

® http:/ /supernova. lbl.gov/~nugent/nugent_templates. html 
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Fig. 8.— The distribution of relative for the SNIa, SNIb/c and 
SNII samples when analyzed in all cases with la SED based light 
curves. The distribution of SNII data are significantly than the 
same for SNIa. This allows Type II SNe that are contaminating 
la samples to be easily distinguished. However, the peak of the x 
distribution for the SNIb/c overlaps with the same for SNIa and 
we find SNelb/c to be the main contaminants for photometrically 
selected la samples for cosmological measurements. 

We find that when fitted with la SED light curves, the 
values for SNell are so large large values that they 
are easily distinguished from the SNIa even with photo- 
metric data alone. However, for SNelb/c the peak 
overlaps with that of the SNela, so they are the primary 
contamination to the total sample. 

To obtain a less-contaminated sample, as a first cut 
we note that the distribution of the SNelb/c has 
a long tail which can extend to tens of thousands, 
and some of the SNelb/c can be removed by an over- 
all restriction on the x^ values in the fitting to SNIa 
light curve template. If the selection is restricted to 
xlei < 20, keeping all real Type la's, this results in a 
removal of about 40% of the SNelb/c's. Since the ra- 
tio of th e rate of SNIa to SNIb/c out to z ~ 1 is about 
10 to 7 (iCalura fc Matteucc i' '2006: 'dc Plaa ct al.' 20071 
Sato et al.ll2007t iPoz nanski et al. .2007a: Eldridge et aTl 
20081; LSmartt ct al. 2009: G eorgv et al.ll2009( ). we expect 
about 250 SNelb/c to remain and contaminate a sam- 
ple that contains 1000 SNela selected photometrically. 
Note that here we just use a simplified assumption with 
the ratio of the rate of SNIa to SNIb/c to be redshift 
independent. In the next subsection we discuss a sta- 
tistical method to further reduce the contamination of 
Ib/c's during model fits to the total sample. 

4.2. The Bayesian Statistical Method 

We employ t he B ayesian estimation method proposed 
bv lPressI (|1996f ) and lKunz et"all (|2007[ ) to further reduce 
the contamination from SNIb/c. We note that our pro- 
posed statistical method cannot distinguish each SNIb/c 
from a Type la individually, but statistically it reduces 
the overall contamination and the associated bias in cos- 
mological parameters. As we illustrate here, the same 
method also allows us to jointly estimate the fraction of 
the SNIa (or Ib/c's) within the whole supernovae sample 
used for cosmology. 
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We take the case that the observational sample of sup- 
posedly Type la's D contains a mixture of true SNela 
data Di and SNelb/c's Dq which mimicks la. We define 
a vector v of length the total number of SNe N with the 
value of Vi taking either 1 or if Di is or is not a SNIa. 
We also define a quantity p to account for the total frac- 
tion of the true SNela in the total SNe data set D. Using 
V and p, the posterior probability can be written as 

P(0|D)=^P(0,v,p|D) (10) 
cx^/:(D|0,v,p)F((?,v,p) (11) 

p,v 

^Y.C{-D\0,v,p)P{p)P{0\p)Pi^\e,p), (12) 



In Eq. ([TOl) , the sum over p will be the integration if the 
value of p is continuous, and the sum of v goes through 
all 2^ possible values of v. Eq. (fTTj) is derived from the 
Bayes theorem, and £(D|6',v,p) is the likelihood. The 
P{0\p) in Eq. p2|) can be reduce to P{0) since there is 
no reason to believe the parameter p affect the cosmo- 
logical evolution of the Universe (it is not a cosmological 
parameter). 
Thus, we can simplify to 

P{0\B) (X P{0) J2 Pip) E '^(I^l^' v,p)/'(v|^,p). (13) 

p V 

For any value of p, P(v|0, p) is when Vi involving the ith 
datum is a SNIb/c. When normalized, P{vi = l\0,p) = p 
and P{vi ~ O\0,p) = 1 — p. Therefore, we find 



p{0\B) (X p{0) p{p) y [ n cIp Y[c°{i~ p) 



p 



(14) 

Here C] is the likelihood that the ith SN is a la and this 
is taken to be 

1 



2'Kai 



(15) 



where ai is the error and Xi = [l^i ~ IJ^ifl^i-, and 
the pf"^ and are the observational and theoretical 
distance modulus respectively. 

The L\ is the likelihood for the SNIb/c samples though 
we don't know a priori the exact distribution. We use two 
parameters h and ctq belonging to the parameter set to 
describe as 



1 



(16) 



/27rao 

where = (/if ^ - - ^Yl<yl- 

We can simplify Eq. (fT4| further by noting that the 
1^ summation term can be written as the product of N 
terms. We finally get 

p[0\D) (X p(0) y Pip) n [cIp + - P)] . (17) 

p N 

The sum over p is easily performed with MCMC runs, 
and we assume P{p) is a uniform distribution. When 
analyzing our mock samples we take the ranges for p, b 
and CTo ofpe (0.5,1), be (-20,20), I/ctq G (0,1000). 



4.3. The Results 
Extending the discussion in §4.1, we add 250 SNIb/c 
data with xlei < ^0 to the 1000 SNIa data, and extract 
cosmological constraints on the time-evolving equation of 
state of dark energy with an analysis which implements 
the Bayesian estimation method described above, in ad- 
dition to a method where all data are analyzed with a 
MCMC run without making an attempt to account for 
Ib/c contamination to the total sample. 
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Fig. 9. — The contour maps for wq vs. wi with and without 
Bayesain estimation. The la, 2cr and 3cr errors are shown. 




Fig. 10. — The PDF of the fraction p for the SNela in the mixed 
sample of la's and Ib/c's. The fiducial value is 0.8 while the fitting 
result based on the technique outlined in § 4.2 is p = 0.75j 



-+0.05 
-0.06- 



We show the contour maps of Wq vs. Wi in FiglHl The 
red solid and blue dotted contours are the fitting results 
with and without Bayesian estimation, respectively. As 
can be seen in Figl9l the constraint on wq and wi is com- 
pletely wrong when we ignore the Ib/c contamination in 
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the total sample and just do direct fitting to the Hubble 
diagram. This is caused by the large differences of the 
distribution between the SNIa and the SNIb/c data and 
MCMC chains are easily trapped in a wrong likelihood 
value. The value for the overall best-fit in this case is 
also very large reaching as high as 2000. 

When we implement the Bayesian estimation method, 
the result is improved significantly. Although there is a 
difference between the best fit and the actual (fiducial) 
value, the fiducial value of the cosmological parameter 
set (wo,wi) = (—1,0) lies safely within the 2ct contour 
around the best fit. Also, as discussed, we also jointly 
estimate the fraction of SNela in the total data set. We 
plot the likelihood for P{p) in FigfTOl The fraction of the 
SNIa in this particular mock data should be 80% while 
the fitting leads to the result of p = O.TSlg Qg with errors 
at la. While there still remains a bias associated with 
the contaminating Type Ib/c's, we have reduced this bias 
to the level of a few percent. Also, if the distributions of 
b and (Jq are better measured in the future, this method 
would get better typing result. 

We believe the Bayesian estimation method provides 
a useful statistical tool to reduce contamination and to 
evaluate the fraction of the contaminating supernovae 
in the LSST photometric SNIa survey. Of course, to 
identify whether each individual SN is a core-collapse 
one or a la, the method discussed above is inadequate, 
but the Bayesian statistical analysis can be a valu- 
able g uide for further advanc e d studv (iPoznanski et al.l 
200l IGal-Yam et al. I [200l IJohnson fc CrottsI 12001 
Poznanski et al. I l2007bt iKlITnetsova fc Connollvll2007n 



5. SUMMARY 

In this paper, we explore the ability to determine the 
redshift and the other parameters useful to construct the 
Hubble diagram with light curve for the LSST SNIa pho- 
tometric measurements. Using a SNIa SED template and 
the expected photometric error of the LSST, we first 
simulate the observed flux data of 1000 SNIa in each 
of 5 LSST filters, and then apply a MCMC technique 
to fit the redshift, stretch factor, apparent magnitude 
and the phase of the SNIa, among others. We find that 
when z < 0.2, these parameters can be determined ac- 
curately at a level comparable to the case where spec- 
troscopic redshift is known. At higher redshifts, the un- 
certainty in photometric redshift goes up quickly since 
o'phot ^ lO'^'^'", but the photometric data is still very 
useful when 0.2 < z < 0.5. To illustrate the effect of the 
uncertainty of the photometric redshift on the fitting of 
the cosmological parameters, we also extract cosmologi- 
cal constraints using parameters of the SNIa light curves 



with and without spectroscopic redshifts. Using the fit- 
ting results of the two cases, we constrain the cosmologi- 
cal parameters for rimo and Q\o in the ACDM model and 
wq and wi in the time-evolved wCDM model. We find 
that for the same number of SNIa data, the cosmology 
fitting with only the photometric data leads to a factor 
of 2 degradation in error of cosmological parameters or a 
factor of 4 in the figure of merit of dark energy equation 
of state (i.e. the inverse area of the Wq — Wi ellipse) com- 
pared with the case of fitting with spectroscopic data. 
However, as the number of photometric-only data far ex- 
ceeds that with spectroscopic data, the overall statistical 
uncertainty in the former would still be smaller. 

Finally, we discuss the contamination on the SNIa 
data from core-collapse supernovae involving types II and 
Ib/c, and the feasibility of using a Bayesian estimation 
statistical method to reduce the overall contamination. 
Similar to SNIa mock samples, we generate the mock 
flux data for the SNIb/c and SNII based on their spec- 
tral templates, and use the SNIa fitting process to fit 
them. We find that the SNell are easily distinguished 
from SNIa because there is an apparent mismatch (large 
X^) when fitting with the SNIa templates. However, this 
is not the case for Type Ib/c's. The peak of its distri- 
bution is overlapping with that of the SNIa and present 
a significant contamination of any photometric selected 
supposedly SNIa samples, even if a conservative cut is 
applied in the values for selection. To further account 
for this contamination, at least statistically when doing 
cosmological model fits, we employ Bayes theorem. Our 
suggested method could reduce the contamination down 
to a few percent level, leading to estimates of cosmolog- 
ical parameters that are biased within la errors. The 
method also establishes the fraction true SNIa in the 
total photometric SNe data set. Nevertheless, we must 
note that this method cannot distinguish if an individual 
SN is whether Type la or not. We will need an extended 
analysis complemented with additional observations if we 
are required to recognize the type of individual SNe. 
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